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The  investigations  reported  in  the  Third  Semi-Annual  Technical  Report 
(Report  No.  16)  have  been  continued  and  extended.  Belcw,  we  describe  the 
research  done,  according  to  the  items  listed  in  the  contracts 

1.  Strain  release  studies. 

2.  Laboratory  studies  of  scale  models  under  stress. 

3*  Development  of  a  method  to  measure  stress  variations  in  the  crust.  < 

In  addition,  the  report  contains  the  following  two  items! 

4.  S.J.  Duda's  visit  to  USA  for  study  and  research. 

5.  Bibliography. 

Complete-  .uscripts  of  two  papers  are  attached  as  Appendixes  I  and  II. 
All  the  research  within  this  contract  has  been  conducted  by  the 
following  two  personal 
Seweryn  J.  Duda,  full-time 


Markus  BAth,  part-time. 
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1 «  Strain  release  studies 

a.  Western  part  of  the  North  American  Continent 

Our  strain  release  studies  are  concerned  with  the  ciroum-Pacific 
earthquake  belt*  The  purpose  is  mainly  to  investigate  possible  mutual 
dependences  of  seismic  activity  in  ad,'  i.oent  areas*  Earlier ,  three  sequences 
in  California  (Kern  County  1952,  Desert  Hot  Springs  1 94S *  San  Francisco 
1957)  and  the  Chile  earthquakes  I960  were  investigated.  See  1*  and  2*  in  the 
Bibliography.  In  this  period,  the  western  part  of  the  North  Amerioan  Continent 
was  included  in  these  investigations. 

First,  the  seismicity  of  West  Canada  for  the  interval  1955-195 9  was 
investigated.  In  the  formula  log  N  »  a  -  bM,  where  N  is  tl  _  number  of  shocks 
in  the  magnitude  range  M  +  0.05 »  we  found  that  the  coefficient  b  -  O.58  for 
M  >2.4*  his  is  a  relatively  low  value  of  b.  One  of  the  main  difficulties 
in  investigations  of  this  type  is  to  limit  off  the  observational  material 
correctly  in  space  and  time.  In  this  case,  we  gradually  found  it  necessary 
to  extend  the  observations  to  cover  the  whole  of  western  North  America 
between  latitudes  10°N  and  67&°N  and  to  extend  them  in  time  to  cover  the 
11-year  interval  1949-1959*  Our  listed  material  is  most  probably  complete 
for  earthquakes  with  M  >6.0.  A  time-space  diagram  has  been  constructed, 
like  the  one  earlier  made  for  Chile  (see  2.).  We  are  still  working  on  the 
search  for  possible  migration  patterns.  Mere  elaborate  statistical  technique 
will  be  used  in  order  to  ascertain  the  significance  of  any  migration  patterns 
found.  Although  there  are  some  indications  of  migration,  there  is  so  far  no 
reliable  proof  of  their  existence  in  this  case. 

b.  Non-AlPide  Asia 

We  have  extended  our  strain  release  studies  to  include  also  one  or  two 
aftershock  sequences  from  the  continental,  non-Alpide  Asia,  mainly  to  get 
observational  material  from  an  area  which,  according  to  some  indications, 
exhibits  a  pattern  different  from  the  Pacific  belt.  We  have  already  begun 
with  the  Mongolian  sequence,  starting  on  December  4,  1957,  and  have 
tentatively  planned  to  include  a  sequence  in  China  (Chinghai  Province, 
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May  21,  1962). 

The  investigation  of  the  Mongolian  sequenoe  from  1937  has  so  far 
yielded  the  following  results.  All  earthquakes  with  magnitude  M  >5.2  have 
been  listed  for  3&  years  following  the  main  shock.  The  b-coefficient  in  the 
formula  log  N  -  a  -  bM  was  found  to  be  o.7  for  this  material,  i.e.  about  the 
same  as  for  Chile  (see  2.).  A  strain  release  characteristic  has  been 
constructed,  which  exhibits  two  branches.  These  can  be  approximated  by  a 
logarithmic  and  an  exponential  ourve  resp.  However,  these  curves  will  be 
reconstructed  in  the  light  of  new  findings  (see  l.d.  below).  In  the  first 
stage  of  the  sequence,  the  epicenters  of  the  aftershocks  were  distributed 
along  gui  E-W  striking  fault.  On  December  3,  I960,  some  seismic  activity 
started  on  an  adjacent  fault,  extending  toward  SE  from  the  eastern  end  of 
the  old  fault,  under  an  angle  of  40°  with  this  fault.  The  activity  on  the 
second  fault  is  most  likely  caused  by  the  pre-existing  activity  on  the  first 
fault.  Our  photoelastic  experiments  have  indeed  provided  an  explanation 
for  this  behaviour  (see  2.  below). 

In  order  to  study  possible  changes  of  focal  mechanism  during  the 
aftershock  sequence,  as  the  case  wsus  with  Kern  Couniy  1952,  reoords  have 
been  collected  from  rnemy  stations  for  a  few  of  the  aftershocks. 

c.  Other  aftershock  sequences 

There  are  quite  a  number  of  other  aftershock  sequences,  which  deserve  to 
be  investigated.  Among  these  might  be  mentioned  two  Kurile  Islands  sequences, 
starting  on  Nov.  6,  1958,  and  Oct.  13,  1963,  respectively.  Such  sequences  will 
be  included  in  the  work  under  the  present  contract  to  the  extent  that  time 
will  permit  and  if  they  are  likely  to  throw  further  light  on  our  problems. 

d.  Improvement  of  the  methods 

Benioff  initiated  strain  release  studies  around  1950.  Since  that  time 
the  same  methods,  as  originally  given  by  Benioff,  have  been  used  by  all  who 
have  worked  in  this  field,  including  ourselves.  In  the  original  method, 
strain  is  proportional  to  the  square  root  of  the  released  energy.  The  volume 
of  every  aftershock  was  considered  constant  and  equal  to  the  total  volume 
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of  the  aftershock  zone.  Moreover,  an  older  energy-magnitude  formula  has 
been  used  for  consistency  reasons,  although  newer  and  better  formulas  have 
been  developed  in  the  meantime.  The  fraction  of  elastic  strain  energy 
converted  into  seismic  energy  was  also  assumed  constant.  It  occurred  to  us 
that  time  had  come  to  improve  Benioff's  method  in  the  light  of  new  findings. 

Under  the  assumption  that  the  volume  of  an  earthquake  is  defined  by  the 
total  volume  occupied  by  its  aftershocks ,  it  was  possible  to  deduoe  a 
relation  between  earthquake  volume  and  magnitude.  The  volume  was  found  to 
increase  with  magnitude,  at  about  the  same  rate  as  energy  increases  with 
magnitude.  Strain  was  found  to  be  independent  of  magnitude,  within  our  error 
limits.  Therefore,  the  difference  between  large  and  small  shocks  is  mainly  a 
difference  in  the  ''olume  and  not  in  the  strain,  a  result  which  may  seem 
surprising  but  is  in  better  accord  with  general  inferences  from  rock  behaviour 
under  stress.  As  a  oonsequence,  some  of  our  earlier  strain  release 
characteristics  have  been  reconstructed,  now  as  deformation  characteristics 
instead. 

A  paper  has  been  prepared  for  publication,  and  ics  complete  manusoript 
is  attached  to  the  present  eport  (Appendix  I). 

2.  Laboratory  studies  of  scale  models  under  stress 

Our  laboratory  studies  aim  at  a  duplication  of  natural  stress  phenomena 
by  means  of  controlled  experiments  in  which  various  parameters  can  be  varied 
at  will  (both  under  static  and,  if  possible,  also  under  dynamic  conditions). 

In  the  autumn  of  1P62,  it  occurred  to  us  that  the  photoelastic  method  would 
be  most  suitable  for  this  kind  of  work. 

Photoelastic  measurements  require  long  experience  in  order  to  lead  to 
meaningful  results  and  a  well-equipped  photoelastic  laboratory  also  means 
quite  a  lot  of  expenses.  However,  we  were  very  fortunate  in  obtaining  the 
expert  advice  of  Dr.  R.  Hiltscher,  Stockholm,  who  is  a  well-known  authority 
on  photoelasticity,  also  we  had  the  fortunate  possibility  to  work  in  the 
laboratory  of  the  Royal  Water  Power  Board  ("Statens  Vattenfallsverk")  in 
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To  date,  the  two-dimensional  stress  field  har,  been  determined  in 
different  plates <  first  in  one  with  an  open  slit,  then  in  a  plate  with  a 
closed  slit  (zone  of  weakness)  under  t.*o  different  angles  to  the  direction 
of  the  applied  pressure  field.  Part  of  the  evaluation  of  the  measurements 
was  made  on  a  digital  computer  of  the  Swedish  Board  for  Computing  Machinery, 
Stockholm.  The  preliminary  investigations,  so  far  carried  out,  are  reported 
with  all  details  in  an  attached  paper  (Appendix  Xl).  Some  features  of 
earthquake  occurrences  are  explained  in  the  light  of  our  new  findings  from 
these  experiments. 

The  steady  stress  state  in  a  perfect  elastic  model  can  only  give  a 
partial  analogy  to  real  tectonic  conditions,  as  the  stress  state  in  the 
orogenetic  shell  changes  with  time  and  the  material  in  the  orogenetic  shell 
exhibits  elastic  afterworking.  Later,  we  will  try  to  construct  models  with 
time-dependent  elastic  properties  and  to  investigate  them  by  the  photoelastic 
analogue  processing. 

5.  Development  of  a  method  to  measure  stress  variations  in  the  crust 

Direct  measurements  of  stress  in  natural  rock  under  a  certai*  time 
interval  in  order  to  follow  its  variations  with  time  are  of  great  importance 
in  connection  with  our  studies.  To  our  knowledge,  the  sond  developed  by 
Prof.  N.  Hast,  Stockholm,  provides  one  of  the  best  possibilities  for  this 
kind  of  research.  Prof.  Hast  originally  developed  his  sond  for  stress 
measurements  in  concrete,  but  has  in  recent  years  applied  it  mainly  to  test 
rock  stability  in  mines.  He  has  collected  a  number  of  instantaneous  point 
measurements,  especially  in  the  Fennoscandian  area.  His  results  are 
certainly  of  great  geophysical  interest. 

The  Fennoscandian  area  is  very  quiet  from  the  seismic  point  of  view, 
and  what  we  could  find  from  continuous  stress  measurements  could  be  stress 
variations  due  to  earth  tides  and  to  3cme  other  phenomena  (variations  of 
atmospheric  pressure  etc).  Such  measurements  should  be  collocated  with 
similar  measurements  in  seismically  active  areas,  especially  in  Japan, 
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where  they  would  be  of  great  importance  also  for  the  earthquake  predicition 
program  which  is  under  way  there.  Prof.  Hast  has  himself  also  tried  to  get 
UNESCO  interested  in  supporting  such  measurements. 

We  have  repeatedly  been  in  contact  with  Prof.  Hast,  who  always  shows 
great  interest  for  the  geophysioal  applications  of  his  method,  but  so  far 
we  have  unfortunately  not  been  able  to  activate  him  sufficiently  to  get 
some  real  collaboration  started.  At  present,  it  seems  to  be  questionable 
to  what  extent  this  will  succeed  during  the  remaining  18  months  of  the 
contract. 


4.  S.J.  Duda*s  visit  to  USA  for  study  and  research 

As  originally  planned,  Ur.  S.J.  Duda  will  spend  some  time  in  USA  to 
oontinue  his  study  and  research  within  the  items  of  the  contract.  We  have 
now  planned  that  Duda  will  spend  the  first  six  months  of  1964  in  USA, 
mainly  at  the  Seibrnological  Laboratory,  Pasadena,  California,  where  he  will 
oontinue  his  researoh  under  the  guidance  of  Prof.  H.  Benioff.  As  already 
agreed,  Duda  will  receive  the  travel  expenses  from  the  contract,  and  the 
cheapest  possible  air  transportation  will  be  used.  As  a  partial  compensation 
for  the  higher  living  costs  in  USA  as  compared  to  Sweden,  I  suggest  that 
Duda's  monthly  payment  be  raised  to  2500*-  Sw.  Cr.  during  his  stay  in  USA, 
all  wiuhin  the  contract. 

5.  Bibliography 

This  bibliography  gives  a  complete  list  of  papers  produced  up  to  now 
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INTRODUCTION. 

Benioff  initiated  strain  release  studies  around  1950.  Since  that 
tine  the  same  method  as  originally  given  by  Bonioff  (1951 )*  has  been 
used  by  all  who  have  worked  in  this  field,  including  ourselves.  In  the 
original  method,  3train  is  proportional  to  the  square  root  of  the 
released  seismic  energy.  The  volume  of  every  aftershock  was  considered 
constant  and  equal  to  the  total  volume  of  the  aftershock  zone.  The 
fraction  of  elastic  energy  converted  into  seismic  energy  was  also 
assumed  constant.  Moreover,  an  older  energy-magnitude  formula  has 
been  used  for  consistency  reasons,  although  newer  and  better  formulas 
have  beer,  developed  in  the  meantime.  In  the  present  paper  an  effort 

I 

is  described  to  improve  Benioff* s  method,  especially  in  the  directions 
mentioned. 

i 


NOTATION. 

We  shall  be  using  the  following  notation  throughout  the  present 
papers 


D 

E 

E1 

e 

F 

H 

h 


=  constants,  cm; 

a  dip  angle  of  the  fault  plane,  degrees; 

=  deformation,  cm"*; 

®  seismic  wave  energy,  ergs; 

■  seismic  wave  energy  of  the  largest  aftershock  in  a 
bequence,  ergs; 

■  average  strain  in  the  focal  region  of  an  earthquake; 

2 

*  fault  plane  area,  cm  ; 


T3 


vertical  extent  of  aftershock  zone,  cm; 
vertical  extent  of  fault  plane,  cm; 
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log 
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£  V 
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**  width  of  fault  plane i  cm; 

«  elastic  strain  energy,  ergs; 

2 

=  r'gidity,  dynes/cm  ; 

=>  length  of  aftershock  zone,  cm; 

■  length  of  fault  plane,  cm; 

■  logarithm  to  the  base  10; 

■  earthquake  magnitude,  equivalent  to  Mg  in  Gutenberg  & 
Richter's  (1956)  notation; 

■  magnitude  of  largest  aftershock  in  a  sequence; 

■  seismic  gain  ratio  *  E/j  (same  as  "loss  ratio"  of 
Lomnitz  1963) J 

2 

=  aftershock  area,  cm  ; 

*  time  interval  after  the  main  shock,  days; 

=  earthquake  volume,  identified  with  total  aftershock 
volume,  cm^j 

■  volume  of  major  part  of  elastic  strain  energy  content,  cm^; 

=  V  -  V  ; 
e 

■  volume  of  tectonic  stress  field  (Lomnitz  1965)*  cm''; 

=  v/V; 

■  width  of  aftershock  area,  cm. 


EARTHQUAKE  VOLUME  AND  MAGNITUDE. 

The  earthquake  vol'jne  is  not  accessible  for  direct  measurements 
and  it  is  natural  that  slightly  different  definitions  have  appeared  in 
the  literature.  Bullen  (1953*  1955)  identifies  the  earthquake  volume 
with  the  strained  region,  in  which  the  material  is  near  breaking-point 
prior  to  an  earthquake.  Bullen  (1965)  uses  the  term  focal  region  to 
denote  the  volume  from  which  the  major  part  of  the  energy  is  issued. 
Tsuboi  (1956)  also  defines  the  earthquake  volume  as  the  one  where  the 
seismic  energy  was  stored  before  the  earthquake.  Gzovsky  (1962)  calls 
"the  spaco  around  the  fracture  in  which  a  redistribution  of  elastic 
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deformation  energy  is  taking  place"  the  earthquake  focus.  Benioff  ( i 955 ) 
defined  the  original  strain  zone  as  the  total  volume  of  the  aftershocks. 
Benioff  (1962)  estimated  on  the  basis  of  some  results  by  Byerly  and 
DeNoyer  ( 1 95Q )  that  the  strain  is  confined  to  a  very  narrow  zone  around 
the  fault,  However,  inferences  from  geodetic  measurements  may  be  of 
limited  applicability  because  of  their  necessary  limitation  to  the  earth's 
surface.  Unfortunately,  there  are  no  measurements  available  which  permit 
a  collocation  of  -;he  ideas  of  t.ie  two  Benioff  papers  mentioned. 

We  define  earthquake  volume  as  the  volume  of  major  energy  content 
and  assume  this  to  be  identical  with  Benioff 's  (1955)  original  strain 
zone.  Table  I  summarize^  all  pertinent  information  on  aftershock 
sequences  we  could  find  in  the  literature,  partly  revised,  and  Pig.  1 
shows  the  corresponding  plot  of  log  V  versus  M.  Thsre  is  no  doubt  that 
the  earthquake  volume  increases  with  magnitude  and  the  straight  line 
in  Fig.  1,  corresponding  to  the  "'ollowing  least-square  solution, 
represents  the  data  wells 

log  V  =  (9,58  +  0.51)  +  (1.47  +  0.14)  M  (1) 

for  5.3  -  M  t  6.7 

It  is  interesting  to  note  that  V  increases  with  M  at  about  the  same 
rate  as  the  seismic  energy  E  does  (B&th  1958)*  We  will  assume  eq.  (l) 
to  be  of  general  validity,  e.g.  also  for  the  individual  shocks  in  an 
aftershock  sequence. 

The  scatter  in  Fig.  1,  reflected  in  the  mean  errors  in  eq.  (l), 
has  several  reasons,  which  we  summarize  as  follows* 

1 .  Casual  errors  of  M  are  likely  to  be  less  than  1/4  of  a  magnitude 
unit,  and  only  of  small  consequence  for  our  relation. 

2.  Errors  of  7  may  be  of  greater  consequence,  because  of  uncertainties 
in  length,  width  and  depth  of  aftershock  zones.  Howe’-  v,  as  we  are 
naturally  only  concerned  with  orders  of  magnitude  of  the  volume,  these 
uncertainties  are  also  only  of  minor  importance.  In  the  sequences 
used,  there  are  some  variations  in  tha  magnitude  range  and  time  interval 
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conaidered,  but  these  have  practically  no  influence  on  the  result. 

Recent  observations  of  free  oscillations  of  the  earth  and  of  surface  waves 
have  permitted  an  estimate  of  the  fault  length  (Benioff,  Press  and 
Smith  1961 )  Press,  Ben-Menahem  and  Tokscz  1/6l,  Bcn-Menahem  and  Toksbz 
1962),  indicating  that  the  length  of  the  aftershock  area  exceeds  the 
original  fault  length  only  by  some  10  percent,  agreeing  with  independent 
data  for  Kern  County  1952,  a  systematic  error  of  no  consequence  here. 
Information  from  geological  expeditions  is  also  in  good  accord  with 
our  assumption.  The  percentage  error  of  the  vertical  dimension  is 
somewhat  greater  than  for  the  other  two. 

3.  Systematic  variations  of  V  may  exist  in  comparing  different 
earthquake  regions.  Unfortunately,  very  little  information  is  available 
on  this  point,  except  for  a  few  hints.  Duda  ( 1 963}  found  indication 
for  a  more  brittle  rheological  behaviour  in  Chile  than  in  the  Aleutian 
Islands  or  Kamchatka.  The  comparatively  low  volume  found  for  the  Kern 
County  earthquake  may  be  explained  by  a  reduced  shear  -*trength  because 
of  many  fractures  and  minor  faults  in  the  area  (Benioff  1955)  or  more 
specifically  by  the  intervention  of  the  Edison  fault  into  the  activity 
on  the  White  Wolf  fault  (Duda  and  B&th  1963)* 

It  must  be  understood  that  relation  (l)  is  a  simplification  of 
real  conditions,  at  least  to  the  same  extent  that  any  relation 
between  E  and  M  is  a  simplification.  A3  indicated  in  point  3*  above, 
there  are  certainly  a  number  of  other  factors  entering  any  such 
relation  as  well,  but  at  present  these  are  impossible  tn  take  into 
account . 

Moreover,  it  must  be  emphasized  that  even  in  case  the  aftershock 
total  volume  V  would  be  essentially  larger  than  the  volume  where 
the  major  part  of  the  elastic  strain  energy  is  stored,  i.e. 

v  .  Ve  +  £V  (2) 

cur  relation  is  correct  up  to  a  constant  factor,  provided  £iV/V  is 
independent  of  magnitude. 
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FAULT  PLANE  AREA,  AFTERSHOCK  AREA  AND  MAGNITUDE. 

Before  proceeding  we  shall  also  consider  some  other  relations 
between  geometrical  properties  of  an  earthquake  and  its  magnitude. 

Berckhemer  (1962)  gave  a  relation  »e tween  fault  plane  area  and 
magnitude ; 

log  F  =  0.45  +  1.7  M  (3) 

for  5*5  -  M  £  8.0 

2 

Although  his  data  would  suggest  introduction  of  an  M  -term  in  this 
relation,  especially  for  M<6,  where  unfortunately  the  dtisa  are  very 
scanty,  only  the  linear  relation  (3)  is  given  in  Berokhemer's  paper. 

Th.j  aftershock  area  S  has  been  related  to  the  magnitude  of  the 
main  shock,  especially  by  several  Japanese  seismologists.  The  latest 
published  relation,  krown  to  us,  was  given  by  Utsu  and  Seki  ( 1 955 ) » 
based  on  39  aftershock  sequences  in  and  near  Japan: 

log  S  =  5.99  +  1.02  M  (4) 

According  to  Dtsu  ( 1 961 )  this  relation  is  confirmed  by  aftershock 
sequences  from  other  areas  as  well,  except  for  the  Kamchatka  1952 
sequence,  which  had  an  exceptionally  large  area.  Fig.  1  shows  the 
straight  line  (4)  together  with  some  of  our  observations.  The  agreement 
is  not  very  satisfactory,  and  a  new  relation  was  derived  fr~a  our 
observations  by  the  least-square  method: 

log  S  -  (4.95  +  O.43)  +  (1.21  +  0.18)  M  (5) 

This  relation  is  certainly  based  on  only  six  observations,  which, 
however,  represent  a  number  of  quite  different  earthquake  regions-  and  line 
up  remarkably  well  on  a  straight  line,  again  wi'bh  exception  for  Kern  County 
1952,  which  is  low. 

From  the  geometry  of  Fig.  2,  showii  1  the  fault  plane  of  the  main 
shock  inside  the  total  aftershock  volume,  in  the  shape  of  a  paralleiopiped, 
we  have  immediately  the  following  relations: 

V.LIH  |  S  =  L  W  j  F  =•  £h1  j  h1  »  h/sf  0C  (6) 

For  the  cases  with  available  information  onoC(Table  J.),  it  is  obvious 
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that  since  2  0.87,  and  therefore  we  can  put  approximately  equal  to 
h»  For  obvious  reasons,  F/LH  £L\ .  This  ratio  is  plotted  against  magnitude 
in  Fig.  1  (see  also  Table  I).  There  is  evidently  some  indication  of  an 
increase  of  F/LH  with  M,  approaching  unity  for  the  largest  shocks,  as 
shown  by  the  dotted  curve  in  Fig.  1. 

ELASTIC  STRAIN  ENERGY  AND  SEISMIC  ENERGY. 

Combining  eq.  (l)  with  BAtb’s  (1958)  energy-  magnitude  formula: 
log  E  -  (12.24  +  1.35)  (1.44  +  0.20)  M  (7) 

we  find  that 

log  (E/V)  =  (2.66  :  1.86)  -  (0.03  +  0.34)  H  (8) 

This  means  that  the  seismic  energy  density,  i.e.  the  seismic  energy 
per  unit  volume,  is  independent  of  the  magnitude.  Fig.  1  shows  the 
straight  line  (8)  together  with  our  observations.  Similar  results  have 
been  expressed  by  Pshennikov  ( 1 962 ) • 

From  the  values  in  Table  I  we  find  that  for  circum-Pncifie 
earthquakes,  cases  i-6,  the  average  value  of  log  (e/E^)  a  1.73  +  0.29, 
corresponding  to  an  average  of  M  -  M1  -  1 .2  +  0.2,  an  excellent 
confirmation  of  the  so-called  B&th's  law  (Richter  1958,  P*  69).  This 
agrees  with  No.  10  (Outer  Mongolia). 

The  source  from  which  earthquakes  derive  their  energy,  i.e. 
the  potential  energy  within  th.  earth,  is  not  accessible  to  direct 
measurements.  The  potential  energy  consists  mainly  of  the  elastic 
strain  energy,  even  if  some  other  kinds  of  potential  energy,  due  to  the 
state  of  the  material,  cannot  be  excluded,  especially  at  greater  depth 
(Benioff  1963).  So  far  we  have  had  no  information  on  the  seismic  gain 
ratio  q  =  E/j  and  its  possible  relation  to  magnitude.  Benioff  (1951 ) 
assumed  q  to  be  constant  and  put  it  equal  to  unity. 

Lomnitz  (1963)  expressed  the  hypothesis  that  q  is  related  to  F 
and  v  in  the  following  way: 


(9) 
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with  A'  and  A  =  A'/p  assumed  by  us  to  be  independent  of  the  magnitude. 

Even  if  the  expression  (9)  seems  to  be  reasonable,  we  have  to  remember 

that  it  is  nothing  more  than  a  hypothesis.  Then,  q  wou3d  increase  with 

magnitude  at  the  same  rate  as  P/V  (Table  II).  Assuming  that  q  has 

reached  its  maximum  value,  i.e.  unity,  for  the  largest  shocks  (here 

taken  as  M  =  8.7),  as  inferred  from  our  result  concerning  magnitude 

7 

variation  of  F/LH,  we  get  the  maximum  value  of  A  *  1.35  *  10  cm. 

With  this  value  of  A  together  with  eqs  (9),  (3)  and  (l)  we  find  that 
log  q  =  (-2.00  +  0.51)  +  (0.23  +  0.14)  M  (10) 

We  see  from  Table  II  that  about  seven  times  more  of  the  elastic  strain 
energy  is  converted  into  seismic  energy  for  an  earthquake  of  magnitude 
8.7  as  for  one  of  magnitude  5*0* 

From  the  definition  of  q  and  eq.  (9)  we  get 

log  J  =  log  E  -  log  A  -  log  F  +  ] jg  V  (ll) 

Inserting  eqs  (l),  (3),  (7)  and  A  *  1.35  x  13^  cm,  we  find  that 
log  J  -  (14.24  +  1.86)  +  (1.21  +  0.34)  M  (12) 

Under  the  same  conditions  we  immediately  derive  the  following 
expression  for  the  elastic  strain  energy  density: 

log  (J/V)  =  (4.66  +  1.35)  -  (0.26  +  0.20)  M  (',>) 

There  is  no  significant  variation  of  J/V  with  M,  just  as  the  case  is 
with  E/V,  eq.  (8) . 

STRAIN  RELEASE  AND  DEFORMATION  CHARACTERISTICS. 

With  reference  to  B&th  and  Benioff  (1958)  we  have  the  following 
formula  for  calculation  of  the  average  strain  £: 

£  =  (14) 

The  strain  is  in  this  formula  composed  of  a  distortional  and  a 
dilatational  part,  which  cannot  be  separated  from  each  other.  Alb 
the  volumes  of  distortional  and  dilatational  strain  energy  storage  nave 
to  be  assumed  equal . 
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In  the  light  of  the  results  described  earlier  in  this  paper,  -we 
suggest  the  following  improvements  of  Benioff's  ( 1 951 )  original  method 
in  the  application  of  eq.  ( 1 4 ) s 

1)  V  will  vary  with  M  according  to  (l),  instead  of  being  assumed  constant} 

2)  q  is  assumed  to  vary  with  M  according  to  (9),  as  an  alternative  to 
the  still  plausible  assumption  of  constant  q  =  1  according  to  Benioff 
(1951)} 

3)  E  varies  with  M  according  to  (7),  being  the  most  reliable  energy- 
magnitude  formula  so  far  produced,  in  excellent  agreement  with  an 
independent  recalt  by  Gutenberg  and  Richter  (1956). 

11  2 

The  rigidity  ju.  will  still  be  assumed  constant  =  6  x  10  dyne s /cm  , 

valid  for  the  upper  part  of  the  earth  where  the  earthquakes  considered 
took  place.  However,  there  are  some  indications  that  fX  depends  on  the 
stress  state  of  the  material  (Duda  1962). 

Solving  eq.  (14)  for  £.  and  considering  only  the  positive  root, 
which  may  be  approximately  correct  as  long  as  we  consider  only  one 
earthquake  area  at  a  time,  we  have 


=  l_2S_V 

\q*V  , 


We  consider  two  cases,  depending  upon  the  expression  for  q  chosen. 
F 

l)  q  =  A  ^  ,  i.e.  eq.  (9).  Eq.  (15)  then  becomes 


=  (— — V 
lvip  I 


n 

Applying  eqs  (?)  and  (3)  and  putting  A  =  1.35  x  10  cm,  we  obtain  that 
log  C  =  -(3.41  +  0.68)  -  (0.13  +  0.10)  M  (17) 

2)  q  =  1.  In  this  case,  (15)  becomes 

(,8) 

Inserting  eqs  (l)  and  (7).  we  find  that 

log  £  =  -(4.41  +  0.93)  -  (0.015  +  0.17)  M  (19) 

The  resulting  eqs  (17)  and  (19)  agree  in  the  sense  that  the  strain  has 
no  significant  variation  with  magnitude.  See  Table  II. 
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This  result  may  seem  surprising  at  first  sight.  However,  it  is 
in  better  accord  with  general  inferences  from  rock  behaviour  under 
stress  than  that  strain  should  increase  rapidly  with  magnitude.  The 
essential  difference  between  large  and  small  shocks  is  not  to  be  found 
in  the  strain  release  but  in  the  volume  within  which  a  re? ease  takes 
place  at  the  same  time.  The  rock  can  store  a  certain  amount  of  strain  n 

before  it  breaks.  If  strain  were  magnitude-dependent  this  would  mean 
that  each  rock  should  be  able  to  store  strain  corresponding  to  a 
certain  minimum  earthquake  magnitude.  It  is  then  a  justified  question 
why  in  seismio  areas  not  only  shocks  above  a  certain  magnitude  exist, 
but  also  a  far  greater  number  of  smaller  shocks.  Our  result  of  constant 
strain  but  magnitude-dependent  volume  seems  to  meet  these  problems. 

In  order  to  study  creep  phenomena  in  aftershock  sequences,  Benioff 
(1950  initiated  the  construction  of  strain  release  characteristics. 

Under  the  term  strain,  we  understand  deformation  per  unit  volume,'  by 
virtue  of  the  fact  that  deformation  in  the  neighbourhood  of  any  point 
can  always  be  expressed  as  the  resultant  of  simple  extensions  (called 
principal  extensions)  in  three  mutually  perpendicular  directions  (called 
the  principal  axes  of  strain).  See  Bullen  ^1963,  p«  17)*  We  still 
consider  a  strain  characteristic  as  valuable  in  describing  the 
behaviour  of  the  rock  under  stress,  but  we  are  now  unable  to  construct 
any  such  curve  because  of  the  following  facts: 

1 )  As  strain  is  independent  of  magnitude,  we  should  need  to  know  all 
aftershocks,  especially  the  large  number  of  small  ones,  t°  be  able  to 

construct  a  reliable  curve.  However,  such  information  is  naturally  not 
available . 

2)  Considering  strain,  we  are  only  concerned  with  one  particular  unit 
volume.  It  would  be  incorrect  to  add  strains  from  quite  different 
volumes,  as  would  be  the  case  if  strains  of  many  small  aftershocks  in 
different  parts  of  a  large  aftershock  region  were  added. 

For  these  reasons,  we  have  to  refrain  from  tracing  any  strain 
release  characteristics.  On  the  other  hand,  deformation  characteristics, 
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referring  to  the  whole  aftershock  volume  can  be  traced  because  of  the 
theorem  mentioned  (Bullen  1963,  p.  17)  and  because  neither  of  the 
objections  above  is  applicable  in  this  case.  In  addition,  as  focal 
mechanisms  are  very  similar  or  closely  related  within  one  and  the  same 
earthquake  area  (B&th  1952),  the  principal  axes  of  strain  are 
approximately  conserved  within  such  an  area  and  it  will  be  justified 
to  add  deformations  from  different  parts  of  the  same  area.  The 
deformation  characteristic  gives  a  true  picture  of  the  real  happenings 
in  the  aftershock  zone  and  has  an  obvious  interest  from  the 
tectonophysical  point  of  view.  Of  course,  under  Benioff 's  (1951) 
assumptions  the  strain  and  deformation  characteristics  had  the  same 
shape,  differing  only  by  a  constant  factor. 

Using  eq.  ( 18 ) ,  i.e.  assuming  q  =  1,  we  have 

0  *  *  '  M4  (20) 

which  in  combination  with  eqs  (l)  and  (7)  becomes 

log  D  =  (5.17  +  0.93)  +  (1.46  +  0.1?)  M  (21) 

E,  V  and  D  all  increase  with  M  at  about  the  same  rate. 

Using  eq.  (21 )  on  two  aftershock  sequences  for  which  the  material 
has  been  published  earlier,  i.o.  Aleutian  Islands  1957  (Dud a  1962) 
and  Chile  I960  (Duda  1963),  we  have  constructed  the  deformation 


characteristics  shown  in  Pig.  3. 

The  accumulated  deformation  in  the  aftershock  zones  can  be 
represented  analytically  as  follows: 

Aleutian  Islands  1957: 
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The  2nd-4th  branches  of  the  Aleutian  Islands  deformation 
characteristic  are  only  to  be  understood  as  straight-line  approximations 
for  an  exponential  curve,  extending  over  the  entire  interval  1,36  ^ 
t^1266,  i.e.  corresponding  to  the  second  branch  for  the  Chile  I960 
characteristic.  Our  improved  methods  have  had  the  consequence  that  the 
earlier  published  characteristics  are  somewhat  changed,  and  it  is  very 
interesting  to  note  that  the  behaviour  is  analogous  in  the  two  sequences 
studied  here.  This  behaviour  was  also  evident  for  quite  a  number  of 
aftershock  sequences  studied  by  the  old  method.  Thus,  the  sequences  of 
Long  Beach  1933>  Imperial  Valley  1940  and  Hawke's  Bay  1931  exhibited 
this  behaviour  (Benioff  1951),  also  some  others  reported  in  the  same 
paper  but  interpreted  differently  by  Benioff,  i.e.  Manix  1947,  Nevada 
1932  and  less  clear  Signal  Hill  1933.  Other  examples  are  Kern  County 
1952  (Benioff  1955)  and  San  Francisco  1957  sequences  (Tocher  1959). 

It  remains  to  be  seen  if  the  improved  technique  presented  in  this 
paper  will  bear  out  that  this  is  a  general  behaviour  for  aftershock 
sequences. 

CONCLUSIONS, 

We  can  summarize  the  results  of  the  present  investigation  in  the 
following  points, 

1.  Earthquake  volume,  identified  with  the  total  aftershock  volume, 
increases  with  magnitude  according  to  the  following  equation: 

log  V  =  (9.58  +  0.51)  +  (1.47  +  0.14)  M 

2.  The  ratio  of  fault  plane  area  to  the  vertical  section  through  the 
aftershock  zone,  i.e.  F/LH,  increases  with  magnitude,  approaching 
unity  for  the  largest  shocks. 

3.  The  aftershock  area  increases  with  magnitude  according  to  the 
following  equation: 

log  S  =  (4.95  +  0.43)  +  (1.21  +  0.18)  M 

4.  The  seismic  gain  ratio  is  expressed  as  follows,  adopting  a  suggestion 
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by  Lomnitz  (1963)* 


Under  this  assumption,  q  increases  with  magnitude. 

5*  The  seismic  energy  density,  E/V,  as  vail  as  the  elastic  strain 
energy  density,  J/V,  are  independent  of  magnitude. 

6.  Strain  is  independent  of  magnitude.  Therefore,  the  main  difference 
between  large  and  email  earthquakes  is  not  to  be  found  in  the  strain 
but  in  the  total  volumes  involved,  This  is  in  agreement  with  Tsuboi's 
(1956)  results. 

7*  The  deformation,  i.e.  the  total  strain  in  the  aftershock  zone, 
increases  with  magnitude  according  to  the  following  formula: 

log  D  =  (5.17  +  0.93)  +  (1.46  +  0.17)  M 
i.e.  at  almost  exactly  the  same  rate  as  the  seismic  energy  E  or  the 
volume  V. 

8.  By  means  of  the  improved  method  given  in  this  paper,  some  earlier 
strain  release  characteristics  (Aleutian  Islands  1957  and  Chile  I960 
sequences)  are  reconstructed,  now  as  deformation  characteristics.  It 
appears  likely  that  most  aftershock  sequences  exhibit  similar  deformation¬ 
time  characteristics. 
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SUMMARY 

An  effort  is  made  to  improve  Benioff's  method  for  investigation 
of  strain  release  in  aftershock  sequences.  The  improvements  may  be 
summarized  as  follows: 

1.  Earthquake  volume  increases  with  magnitude,  instead  of  being  constant. 
A  relation  is  given,  relating  volume  to  magnitude. 

2.  A  revised  energy-magnitude  formula  is  used. 

3.  The  seismic  gain  ratio,  i.e.  the  ratio  between  seismic  energy  and 
elastic  3train  enery,  probably  increases  with  magnitude,  instead  of 
being  constant.  Likewise,  the  ratio  of  fault  plane  area  of  the  main 
shock  to  the  vertical  section  through  the  aftershock  volume  increases 
with  magnitude. 

4.  The  seismic  energy  density,  the  elastic  strain  energy  density  as 
well  as  strain  are  independent  of  magnitude. 

5.  The  deformation,  i.e.  the  total  strain  in  the  aftershock  zone, 
increases  with  magnitude  at  the  same  rate  as  seismic  energy  and  volume 
do. 

As  a  consequence  of  these  improvements  some  earlier  published 
strain  release  characteristics  are  reconstructed,  this  time  as 
deformation  characteristics  instead. 
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Table  I  -  DATA  FOR  EARTHQUAKES  USED  IN  THIS  STUDY 
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Table  II  -  VARIATION  OP  GAIN  RATIO,  ENERGY,  STRAIN  AND  DEFORMATION  WITH  MAGNITUDE 
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(Captions  for  the  figures) 

Pig.  1  -  Various  earthquake  parameters  plotted  against  magnitude. 

The  circles  are  our  observations,  using  reference  numbers  of  Table  I. 

Pig.  2  -  Schematical  picture  of  fault  plane  inside  aftershock  lume. 

Pig.  3  -  Deformation  characteristics  for  the  Aleutian  Islands  1957 
and  Chile  I960  aftershock  sequences. 
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THE  STRESS  FIELD  AROUND  A  FAULT 


ACCORDING  TO  A  PHOTOELASTIC  MODEL  EXPERIMENT 
Seweryn  J«  Duda 
(Received  1964  January  ) 

Summary 

A  report  is  given  of  preliminary  photoelastie  measurements  of  the  two- 
dimensional  stress  field  around  a  fault  or  crack  in  a  plate.  The  measurements 
include  the  following  cases:  open  slit  with  the  uniaxial  applied  pressure 
field  making  an  angle  of  45°  with  the  slit;  closed  slit  (zone  of  weakness) 
with  two  different  thicknesses  of  the  weak  zone  and  again  45°  to  the  external 
pressure  field;  the  measurements  for  the  weakest  slit  were  made  also  with  an 
angle  of  22.5°  to  '.he  pressure  field.  The  measured  and  calculated  normal  and 
shear  stresses  are  represented  in  graphical  form.  The  results  provide 
explanations  for  some  earthquake  characteristics,  for  example  distribution 
of  shear  stress  and  patterns  of  geographical  extension  of  seismic  activity 
during  an  aftershock  sequence. 

1.  Introduction 

A  knowledge  of  the  stress  field  around  tectonic  faults  is  of  great 
tectonophysical  importance.  However,  generally  only  indirect  methods  are  at 
our  disposal  for  its  determination.  These  methods  work  with  simplified 
models,  either  mathematical  models  or  laboratory  models.  We  shall  only  be 
considering  two-dimensional  models. 

Inglis  (1913)  gave  the  analytical  solution  for  the  two-dimensional 
stress  state  around  an  elliptical  hole  in  an  infinite  plate  for  different 
'external  stresses.  Anderson  (1951)  applied  these  results  to  a  degenerated 
ellipse,  i.e.  with  one  axis  equal  to  zero,  under  the  following  boundary 
conditions: 

1.  external  pressure  and  tension,  mutually  perpendicular,  act  under  45°  to 
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the  fracture,  so  that  &  pure  shear  stress  arises  parallel  to  the  fracture} 

2,  there  is  no  trr  nsmission  of  shear  stress  across  the  fracture* 

The  principal  stresses  along  the  fracture  were  found  to  be  zero  in  the  middle, 
Increasing  to  infinity  towards  its  ends.  This  is  a  first  approximation  o?  a 
tectonical  fault  in  a  pure  shear  stress  field  and  furnishes  a  good  explanation 
of  some  observed  earthquake  patterns  (Duda  1962),  especially  the  accumulation 
of  epicenters  and  even  more  of  seismic  energy  towards  the  ends  of  a  fault. 
However,  Anderson's  model  is  time- independent,  whereas  an  aftershock  sequence 
continues  for  a  considerable  time.  Therefore,  only  the  time-integral  of  the 
aftershock  pattern  can  be  compared  with  Anderson's  model. 

The  condition  2.  above  is  a  severe  simplification  of  natural  conditions. 

A  real  fault  has  to  be  understood  as  a  zone  of  weakness,  still  able  to  transmit 
shear  stresses  from  one  side  to  the  other.  A  generalization  of  Anderson's 
model  introducing  friction  at  the  fault  would  lead  to  considerable  mathematical 
difficulties,  and  to  our  knowledge  this  problem  has  not  been  analytically 
solved. 

In  the  present  research  we  have  tried  to  solve  this  problem  by  an 
analogue  processing,  using  the  photoelastic  method  for  determination  of  the 
stress  field  around  a  rectilinear,  finite  fracture  in  a  transparent  plate.  Our 
boundary  conditions  are: 

1.  external  homogeneous  pressure  field  acting  parallel  to  the  plate  under 
angles  of  45°  and  22.5°  resp.  to  the  fracture} 

2.  transmission  of  shear  stress  across  the  fracture. 

The  application  of  a  uniaxial,  homogeneous  pressure  field  on  the  plate 
represents  a  more  general  analogy  to  tectonical  conditions  than  a  pure  shear 
stress  field  parallel  to  the  fault.  At  any  point  in  the  orogenetic  shell,  the 
two  principal  stresses  will  in  most  case3  be  two  pressures,  generally  unequal. 
Such  a  pressure  field  may  be  divided  into  a  hydrostatic  pressure  and  a 
residual  pressure,  of  which  only  the  latter  is  of  importance  for  tectonic 
movements.  Its  action  will  depend  on  the  angle  it  forms  with  the  fault.  If 
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the  pressure  is  perpendicular  to  the  fault,  this  will  be  closed  and  no 
earthquake  occurs.  If  they  are  parallel,  the  fault  will  be  opened  or  bent, 
and  no  seismic  activity  will  be  observed  either.  Seismic  activity  will  be 
expected  for  all  other  angles,  being  a  maximum  for  an  angle  of  45°. 

The  photoelastic  method  has  so  far  been  used  only  exceptionally  in 
seismology  or  general  geophysics,  as  e.g.  for  studies  of  seismic  wave 
propagation  (Thomson  1963)  or  of  stress  distribution  in  models  of  wells  and 
galleries  (Schmidt  1961). 


2.  The  photoelastic  method 

Extensive  treatments  of  the  photoelastic  method  have  been  given  by 
Hiltscher  (1958)  and  by  Frocht  (i960).  The  basic  physical  principle  of  the 
method  is  that  some  transparent  materials,  when  loaded,  show  double  refraction 
of  light.  The  velocity  of  light  changes  in  a  loaded  model  from  point  to  point 
and  depends  upon  its  direction  of  propagation.  Thus  a  light  beam  leaving  the 
loaded  model  will  be  polarized.  In  a  two-dimensional  model,  like  a  plate, 
loaded  parallel  to  its  plane,  a  light  beam  will  be  plane-polarized  at  every 
point  in  the  directions  of  the  principal  stresses.  The  two  light  components 
will  have  different  velocities,  proportional  to  the  respective  principal 
stress. 

We  introduce  the  following  notation: 

&  2  ~  PrinciP®l  stresses} 

R  =  relative  retardation  or  phase  difference,  expressed  in  wave  lengths 

C  =  stress-optic  coefficient  (Frocht  I960)} 

d  =  plate  thickness} 

4  d  =  change  of  plate  thickness  due  to  load} 

k  =  lateral  constant  (Hiltscher  1963). 

Then,  the  following  two  equations  are  valid: 
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Rt/d  =  c  (<rx  -<f2)  (i) 

and 

+(f2  =  k4d/d  (2) 

The  elative  retardation  and  the  directions  of  the  principal  stresses 

i 

are  measured  in  a  polariscope  (Figure  l).  Eq,  (l),  called  the  stress-optic 
law  in  two  dimensions,  provides  a  possibility  to  determine  the  difference  of 
the  principal  stresses,  if  R^  is  measured  and  d  and  C  are  known. 

The  sum  of  the  principal  stresses,  expressed  in  eq,  (2),  can  be 
determined  by  different  mechanical,  optical  or  electrical  methods.  In  our 
experiments  we  used  the  very  efficient  and  accurate  lateral  extensometer  after 
Hiltscher  (Figure  2)  for  the  measurement  of  the  thickness  change  Ad.  Knowing 
k  and  d,  we  can  calculate  the  sum  of  the  principal  stresses  from  (2),  Figure  3 
shows  schematically  the  complete  apparatus  used  for  the  photoelastic 
measurements. 

The  described  measurements  and  eqs  (l)  and  (2)  reduce  the  problem  to 
express  stresses  in  any  coordinate  system  to  a  matter  of  calculation. 

3.  Description  of  our  experiments 

Four  different  models,  all  being  transparent  plates  made  of  Araldite  D, 
were  measured  by  the  photoelastic  procedure  (Table  1  rnd  Figure  4).  The  slits 
were  in  every  case  located  in  the  center  of  the  plate  and  orientated  along  a 
diagonal.  Model  0  (where  0  stands  for  "open")  had  an  open  slit,  the  other 
models  close  slits*  consisting  of  cuts  from  either  side,  0.4  cm  deep  in 
Model  A  and  0,45  cr  deep  in  Models  E  and  C.  Profiles  I  snd  II  were  parallel 
to  the  slit  and  profile  V  perpendicular  to  it.  Profiles  III  and  IV,  located 
symmetrically  to  I  aud  II  with  respect  to  the  slit,  do  not  need  to  be  measured 
for  symmetry  reasons. 

Under  load,  the  slit  in  Model  0  is  expected  to  close  and  friction  between 
its  walls  to  arise.  This  corresponds  to  a  fault  after  a  rupture  of  the 
material  along  the  fault  plane.  As  it  is  impossible  to  know,  if  the  slit  is 
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first  closed  and  then  friction  arises  or  if  the  walls  are  first  displaced  and 
then  touch  each  other,  the  three  other  models  were  made.  In  these  the  slit 
corresponds  to  a  zone  of  weakness,  the  strength  being  1/5  of  the  surrounding 
medium  in  Model  A  and  1/10  in  Models  B  and  G. 

In  case  of  Models  A,  B  and  C  there  is  really  a  transmission  of  shear 
stress  across  the  slit.  These  models  correspond  to  pre-existing  but  blocked 
faultB  and  the  stress  field  after  loading  corresponds  to  the  one  preceding 
an  earthquake.  Model  B  gives  the  stress  field  around  the  fault  nearer  the 
yielding  point  as  compared  to  Model  A.  Model  C  is  also  relatively  near 
yielding,,  but  the  shear  stress  is  no  longer  the  maximum  possible  for  the 
applied  pressure,  because  the  angle  between  applied  pressure  and  the  slit  is 
now  only  22.5°  aB  compared  to  45°  in  all  preceding  cases. 


4.  Results 

We  will  be  us Id*  the  following  notation: 
x,  y  =  rectangular  coordinates  in  the  plane  of  the  plate,  the  x-direction 

coinciding  with  the  direction  of  applied  press’ire  (Figure  4)$ 
a(.  -  angle  between  G  ^  and  positive  y-axis,  counted  positive 
counterclockwise  («*<-  9o°)j 


*"o 


=  applied  pressure,  only  absolute  values  used; 

«%■>  tf'v  =  x  8111,1  y  cemponents  of  pressure  (negative)  or  tension  (positive) j 

A  j 


=  reference  shear  Btress,  chosen  =  < T"  * 

’  O' 


fjy  =~Iyx  =  aliear  sto®88:  positive  for  clockwise  shear,  negative  for  anti¬ 
clockwise  shear. 

The  following  formulas  are  us« .  in  the  computation: 


°x  = 


cr,  + 


*1  " 


cos  2  ci 


<rn  +  05 


°y  = 


-  <5? 


cos  2  :/• 


(3) 


T  = 
Lxy 


-  °2 


sin  2 


2 
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The  results  are  given  in  graphical  foru  in  Figures  5-8.  All  stress 
components  are  expressed  in  relation  to  the  absolute  value  of  the  applied 
external  field,  which  renders  all  curves  immediately  comparable.  As  the  slit 
has  a  finite  width  in  Models  A,  B  and  C,  stress  values  near  the  slit  are  less 
reliable  because  of  three-dimensional  effects.  Some  of  the  results  need  to  be 
specially  emphasized  and  the  corresponding  remarks  are  collected  in  what 
follows. 

a.  Profile  V,  Models  0,  A.  B  (Figure  5) 

<j“x  exhibits  a  maximum  value  at  some  distance  from  the  slit,  decreasing 
both  towards  and  away  from  the  slit.  It  is  a  tension,  except  for  Model  0, 
where  it  is  a  pressure  near  the  slit.  The  maximum  value  of  CT^  increases  from 
Model  A  (far  from  yielding)  to  Model  B  (nearer  yielding)  and  is  highest  for 
the  Model  0  (open  slit). 

<3-  is  always  a  pressure  with  a  maximum  value  at  some  distance  from  the 

y 

slit,  decreasing  rapidly  towards  the  slit  and  approaching  cT'  away  from  the 
slit. 

X  has  a  high  positive  value  far  away  from  tne  slit,  and  a  high 

*y 

negative  value  when  approaching  the  slit. 

b.  Profile  I,  Models  0,  A.  B  (Figure  6) 

There  is  no  more  any  symmetry  of  but  the  extremities  exhibit  besides 

double  maxima  a  high  negative  value  to  the  left  in  Figure  6,  and  a  high 
positive  value  to  the  right.  This  stress  distribution  corresponds  to  a 
tendency  of  the  fault  to  rotate  in  the  considered  external  pressure  field. 

This  tendency  is  naturally  most  pronounced  in  the  case  of  an  open  slit. 

Also  cf  exhibits  asymmetry  and  decreases  near  the  middle  of  the  open 

y 

slit  far  below  or^.  There  is  no  axial  symmetry  because  of  the  fact  that  the 
profile  runs  at  some  distance  from  the  slit.  A  parallel  profile  on  the  other 
side  of  the  slit  and  at  the  same  distance  from  it  would  provide  reversed 
mirror- j.mage  curves  to  those  given  here.  This  is  true  in  all  cases  considered. 

The  asymmetry  ofX  with  bigger  absolute  values  to  the  right  in  Figure  6 

*y 


than  to  the  left  for  Models  A,B  is  also  clear  from  Figure  4,  showing  the 
isochromatics  for  Model  A.  Curves  of  equal  colour  connect  points  wi u  equal 
shear  stress,  increasing  from  red  to  green.  A  high  density  of  isochromatics 
corresponds  to  a  rapid  increase  of  the  shear  stress.  In  Figure  4,  Profile  I 
between  points  15  and  21,  there  is  a  region  of  high  density  of  iso chroma tics, 
corresponding  to  the  large  shear  stress  in  Figure  6.  For  the  open  slit,  the 
largest  negative  shear  stress  appears  on  the  left  side  of  tne  profile. 

In  all  models  (0,A,B)  there  are  positive  as  well  as  negative  shear 
stresses  on  the  left  part  of  Profile  I  but  only  negative  stresses  on  the  right 
part.  This  result  may  be  of  importance  in  the  interpretation  of  fault  plane 
solutions  in  aftershock  sequences. 

c.  Profile  II.  Models  0.  A.  B  (Figure  7) 

The  fact  that  Profile  II  is  at  a  greater  distance  from  the  slit  has  the 
result  that  the  stress  curves,  although  similar  to  those  for  Profile  T,  are 
more  smoothed.  For  instance,  the  double  maxima  have  disappeared. 

d.  Profiles  V,  I,  II.  Model  C  (Figure  8) 

Model  C  differs  from  B  only  in  the  angle  between  applied  pressure  and 
the  slit.  Whereas  in  Profile  V,  shows  a  shape  similar  to  Models  0,  A,  B, 
there  is  a  pronounced  difference  in  the  shape  of  the  o'- -curves.  For  the  shear 

y 

stress  along  Profile  V,  the  main  difference  is  the  negative  value  ofT  at 

xy 

the  ends  of  the  profile  in  Model  C  as  compared  to  Models  0,  A,  B. 

The  stress  curves  for  Profile  II  are  more  smooth  as  compared  with  Profile 
I,  just  as  the  case  is  for  the  other  models.  For  these  two  profiles,  all 
stress  values  are  smaller  for  Model  C  than  for  Models  0,  A,  B. 

e .  General  remarks 

As  we  have  measurements  only  along  a  few  profiles,  these  are  naturally 
insufficient  to  give  a  complete  two-dimensional  picture  of  the  stress 
distribution  around  the  fault,  e.g.  in  the  form  of  stress  isolines.  On  the 
other  hand,  the  measured  profiles  are  sufficient  to  give  a  clear  idea  of  the 
main  features  of  the  stress  distribution. 
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We  have  started  from  a  given  external  pressure  field  and  measured  the 
resulting  stress  distribution  around  a  fault.  Applied  to  nature,  the  problem 
may  arise  to  reverse  the  procedure,  i.e.  to  start  from  a  number  of  stress 
measurements  around  a  fault  and  to  deduce  the  external  pressure  field.  This 
is  an  enormously  more  difficult  problem,  and  it  is  not  even  certain  if  it  has 
a  uninue  solution. 

5.  Discussion 

The  purpose  of  this  report  is  l)  to  present  an  efficient  analoque 
processing  method  for  investigation  of  stress  distributions  around  model 
faults  and  2)  to  present  these  distributions  in  a  few  cases  resembling 
tectonical  faults.  The  ultimate  goal  is  to  explain  the  behaviour  of 
earthquakes.  In  this  preliminary  work  we  like  to  point  out  a  few  results, 
which  already  seem  to  provide  explanations  for  some  phenomena  found  in  nature. 

1.  In  the  photoelastic  model  experiments,  with  improved  boundary 
coalitions  as  compared  with  earlier  mathematical  models,  the  increase  of  shear 
stress  towards  the  ends  of  a  fault  is  confirmed.  There  is  also  no  longer  an 
infinite  shear  stress  at  the  fault  ends. 

2.  The  shear  stress  reaches  a  much  higher  absolute  value  at  the  right 
end  of  Profile  I,  Models  A,  B  (Figure  6)  than  on  the  left  end.  For  a  profile 
on  the  opposite  side  of  the  slit,  the  behaviour  is  opposite.  The  strain 
release  density  for  the  Desert  Hot  Springs  aftershocks  (Richter,  Allen  & 
Nordquist  1958,  fig,  6)  exhibits  a  similar  pattern  in  relation  to  the  Mission 
Creek  fault,  which  may  be  explained  by  our  finding.  This  would  mean  that  our 
two-dimensional  stress  distribution  is  a  good  approximation  to  the  stress 
field  around  Mission  Creek  fault,  also  that  the  hypocenters  do  not  vary  much 
in  depth.  Of  course,  the  aftershocks  occur  in  a  material  with  elastic 
afterworking,  and  the  strain  release  density  map  shows  a  summary  effect  ?n 
time,  whereas  the  photoelastic  material  used  by  us  had  to  be  nearly  perfectly 


elastic. 
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3.  In  Figure  4,  the  isochromatics  at  both  ends  of  the  slit  have  a  shape 
resembling  the  wings  of  a  butterfly.  It  is  well  known  from  laboratory 
experiments  that  the  fracture  will  propagate  preferably  in  the  direction  of 
the  "butterfly's  body"  (oral  communication  by  R.  Hiltsoher).  This  direction 
depends  strongly  on  the  external  pressure  field.  In  an  investigation  by  the 
present  author  of  che  Mongolian  aftershock  sequence  starting  on  December  4, 

1957  (03  37  50  GMT,  45.1°N,  99»4°E)  it  was  found  that  the  aftershocks  were 
distributed  along  the  E-W  striking  fault  of  the  main  shock  up  to  a  shock  on 
December  3»  I960  (04  24  18.9  GMT,  42.9°N,  104. 4°E).  With  this  shock  seismic 
activity  started  along  a  fault  extending  in  SE  direction  from  the  eastern  end 
of  the  old  fault,  making  an  angle  of  40°  with  this  fault.  As  explanation  for 
this  secondary  activity  we  propose  that  the  primary  fault  propagated  in  the 
direction  of  the  body  of  the  stress  "butterfly". 
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(Captions  for  the  figures) 

PIG.  1.  -  Polariscope  after  Hiltscher  (1958).  a  =  analyzer,  p  =  polarizer 
with  light  source. 

PIG.  2.  -  Hiltscher  lateral  extensometer  (Hiltscher  1963)  for  measurement  of 
plate  thickness,  pp  =  plate,  mt  =  measuring  tips,  ms  =  micrometer  screw, 
d  =  dial  gauge,  s  =  suspending  device  of  the  extensometer  and  its 
counterweight  c. 

PIG.  3.  -  Loading  and  measuring  device  for  photoelastic  experiments,  after 
Hiltscher  (1958).  Loading  device  to  the  left,  front  view.  Measuring  device 
(polariscope)  to  the  right,  side  view. 

PIG,  4.  -  Isochromatics  for  Model  A  with  closed  slit,  showing  the  profiles 
V,  I,  II  and  the  measuring  points  in  relation  to  the  slit  s. 

PIG.  5.  -  Stress  components  along  Profile  V,  Models  0,  A,  B.  The  measuring 
points  are  indicated  on  the  horizontal  scales,  in  every  case  starting  with 
point  no,  1  at  the  left  end  (compare  Figure  4).  s  =  slit, 

PIG.  6.  -  Stress  components  along  Profile  I,  Models  0,  A,  B.  s  =  slit. 

PIG.  7.  -  Stress  components  along  Profile  II,  Models  0,  A,  B. 


PIG.  8.  -  Stress  components  in  Model  C 


Direction  of  external  pressure 


